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Abstract 

There is a large discrepancy between results of Crystal Barrel and WA102 for the 
branching ratio R = BR[n 2 (1870) -»• a 2 (1320)vr]/ J B J R[?72(1870) -»• / 2 (1270)ry]. An extensive 
' re- analysis of the Crystal Barrel data redetermines branching ratios for decays of 1J2 (1870), 

r/2(1645), r/2(2030) and /2(1910). This re-analysis confirms a small value for R of 1.60=L0.39, 
inconsistent with the value 20.4 ± 6.6 of WA102. The likely origin of the discrepancy is 
that the WA102 data contain a strong /2(1910) — > a 2 ir signal as well as 7/2(1870). There is 
^ ■ strong evidence that the 772(1870) has resonant phase variation. A peak in /2(1270)ao(980) 

CN ■ confirms closely the parameters of the a 2 (2255) resonance observed previously. A peak in 

772(2030)-7r is interpreted naturally in terms of 7^(2245) with reduced errors for mass and 
width M = 2285 ± 20(stat) ± 25(syst) MeV, T = 250 ± 20(stat) ± 25(syst) MeV. 

' i PACS numbers: 11.80.Et, 13.25-k, 14.40. Be 

& 

1 Introduction 

(N ; 

The objective of this paper is to re-examine data from Crystal Barrel (CBAR) [Tj [2] [3] concern- 
ing 772(1645), 772(1870), 7/2(2030) and /2(1910). One motivation is to study a large discrepancy 
for branching ratios of 7/2(1870). Earlier CBAR work found a branching ratio 

g! R= B ^ m \^<i m W=ixr±o.i7, (1) 

O : BR[ V2 (1870) -> / 2 (1270)77] ' 1 } 
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see Table 5 of Ref. j3]. This is much smaller than the value 20.4 ± 6.6 claimed by the WA102 
collaboration in central production of 777^ jl] . Other branching ratios determined by the CBAR 
data are redetermined here. This re-analysis incorporates many further details of spectroscopy 
in this mass range which have appeared since the year 2000. 

A review of these earlier publications will set the scene and introduce the relevant resonances. 
For a comprehensive review of the CBAR work with p in flight see Ref. 0. The first study of 
77 2 (1645) and 77 2 (1870) in CBAR data was at two beam momenta: 1200 and 1940 MeV/c pQ [2]. 
The 772(1645) was observed decaying to ci2(1320)7r and a (980)7r. No other ao or a 2 appear in 
the data, so ao will be used hereafter as a shorthand for a (980) and likewise a 2 for a 2 (1320). 

There was also a strong ^2(1270)77 signal near its threshold ~ 1810 MeV. It could not be ex- 
plained as the high mass tail of 7/2(1645). The reason was that a single 7/2(1645) decaying to f 2 r\ 
would contain a large f 2 r\ signal in both numerator and denominator of the Breit-Wigner ampli- 
tude; cancellation between numerator and denominator cannot accomodate the large f 2 (1270)r] 
signal. The data were fitted with the addition of the 7/ 2 (1870) though the data did not rule out 
the possibility of a non-resonance threshold effect at that time. 
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Figure 1: Scatter plots of M(tt7t) v. M(i]tt) for three ranges of M(r]7r7r): (a) 1560-1750 MeV, 
centred on ^(1645), (b) 1775-1945 MeV over 7? 2 (1870) and (c) 1945-2115 MeV, centred on 
772(2030). The beam momentum is 1350 MeV/c. 



The majority of CBAR data with p in flight were taken in 1996. For 77371"°, statistics were 
a factor 7 higher than earlier data at each of nine beam momenta ranging from 600 to 1940 
MeV/c, i.e. an overall increase of statistics by a factor ~ 30. There are typically 25-30K events 
at each beam momentum and a total of 213K events. 

Fig. 1 shows scatter plots from CBAR data in three ranges of r/Tm mass. Fig. 1(a) shows 
the 772(1645) mass range; there is a vertical band due to g^tt and also a peak in im near 1 GeV. 
The spin analysis ruled out /o(980), showing that the peak in tttt is due to the low mass tail of 
7? 2 (1870) ->■ / 2 (1270)?7. Fig. 1(b) shows the t? 2 (1870) mass range. The a 2 (1320) and / 2 (1270) 
bands appear of similar strength; however, the / 2 is somewhat broader and therefore stronger. 
There is clearly no large excess of a 2 decays in this mass interval. The branching ratio quoted by 
WA102 would require an 7/ 2 (1870) — > a 2 7r signal a factor ~ 3.5 larger than f 2 (1270)7/ when one 
allows for relative decay rates of a 2 — > 7777 and / 2 (1270) — > 7r°7r° and for charge combinations. 
Fig. 1(c) shows the mass range of 77 2 (2030); there is a peak where a 2 (1320) and / 2 (1270) bands 
cross. The a 2 7r band is strong and / 2 (1270)t] is weaker. Note also that weak a (980) bands are 
visible in all three panels. 

There is additional evidence for 77 2 (2030) — > / 2 (1270)t] and a 2 7r in further CBAR data for 
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pp — > r]iT TT , see Fig. 11 of Ref. [5]. A distinctive feature of this state is its strong decay to 
[o27t]l=2; where L is the orbital angular momentum in the decay. 

Following the year 2000 publication of 77371"° results, a combined analysis was made of 10 sets 
of data with isospin 7 = and C-parity +1. Six of these were CBAR data for pp in flight 
— > 7r°7r°, 7]r], rji]', r/7r°7r°, tj'tt ^ and rjrjrj [7]. The other four were high quality differential cross 
sections and polarisations for pp — > 7r + 7r~ from two experiments: Eisenhandler et al. [8] and 
PS172, Hasan et al. [9]. This analysis revealed two towers of resonances centred at ~ 2000 and 
2270 MeV, with all J p for qq states expected in this mass range. Most have been observed in at 
least three sets of data. Polarisation data provide a clean separation of pp 3 P 2 and 3 F 2 states. 
The / 2 (1910) of the PDG [10] was confirmed and identified as dominantly 3 P 2 ; a neighbouring 
3 F 2 state was observed at 2001 MeV. Further / 2 (2240) 3 P 2 and / 2 (2295) 3 F 2 states were also 
observed. In the new analysis reported here, the first three of these f 2 states play a significant 
role. The / 2 (1910) lies close to r] 2 (1870) and is important for the discussion of WA102 data. 

In 2001, a combined analysis was also made of data on 7=1, C = +1 states [TT]. The 
spectrum of states is less complete than for I = 0, C = +1, because of the lack of polarisation 
data. Nonetheless, an a 2 (2255) appeared clearly in three channels of data. Secondly, there is a 
less well identified 7r 2 (2245). These two states now appear in the analysis reported here. So, in 
summary, the picture has developed significantly since the earlier analysis of 77371"°. 

The 7?37r° channel may appear to be a complicated channel to analyse, because of the mul- 
tiplicity of r/7r° and 7r°7r° combinations. However, r/ 2 (1645) and t? 2 (1870) were found here; also, 
once one knows the mass and width of r] 2 (2030) from 7]7nr final states, it is easily detected in the 
present data via its strong decay to [a 2 7r] £=2 , which has a very distinctive angular dependence. 
Interferences between channels provide intricate information on identifiable resonances, even 
though the angular correlations cannot be displayed because they are multi-dimensional. It is 
necessary to work from log likelihood and mass projections of t]tt, tttt, 3tt and rjim. A valuable 
check on the analysis is to introduce deliberate errors into angular dependence of amplitudes; 
genuine signals then drop to low values. 

Two alternative starting points have been adopted. The first is the year 2000 analysis, which 
now improves. The second starts from the WA102 ratio for ?/- 2 (1870) — > a 2 7r and / 2 (1270)r/\ This 
gives a considerably worse fit. After iterations, the two fits converge to a single solution except 
for minor points of ambiguity in a^it decays. No significantly different solutions have been found 
at any beam momentum despite searches from a variety of initial parameters. 

The layout of the paper is as follows. Sections 2 and 3 go through technical details of the 
analysis procedure. It may be useful to glance first at figures of later sections, so as to appreciate 
the rationale for the steps discussed in Sections 2 and 3; the techniques need to be adapted to 
what is found in the data. Section 2 introduces the channels which are required and deals with 
formulae. These are well known from earlier literature, but need to be documented. The one 
point requiring special treatment is the opening of the / 2 (1270)r] threshold, close to r] 2 (1870). It 
is necessary to fold the width of / 2 (1270) into the phase space for the f 2 rj final state appearing 
in the Breit-Wigner amplitude. Secondly, the way the spin dependence is treated in terms of 
partial waves is discussed. 

Section 3 presents features of the data. Fig. 4 shows mass spectra for 37r, tttti], tct] and 7T7t at 
one representative momentum, 1642 MeV/c; other momenta show similar features and Fig. 1 
of Ref. [3] presents spectra at 1800 MeV/c. This section gives more detailed comparisons with 
data in further figures. 
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Section 4 then presents essential results. Table 1 shows changes of log likelihood when each 
channel is dropped from the fit and all others are re-optimised. This Table identifies directly 
the important channels and their significance levels, and how they vary with beam momentum. 
Fig. 7 presents the cross sections for all reactions as a function of beam momentum; these lead 
to considerable insight into the physics. Peaks shown in Fig. 8 may be identified with a 2 (2255) 
and tt 2 (2245). 

A revised set of branching ratios is derived for ?? 2 (1645), 7/2(1870), 7^(2030) and / 2 (1910). 
The strong decay modes change little, but there are some significant changes from the ear- 
lier publication in weak channels, for reasons which are understood. The important ratio for 
7? 2 (1870) — > a 2 (1320)7r and / 2 (1270)r? changes only slightly and remains completely inconsistent 
with the WA102 result. 

Section 5 therefore re-examines WA102 mass projections for r? 2 (1645) and 7? 2 (1870) — > aovr, 
a 2 7r and / 2 (1270)?7. There is good agreement with CBAR data for t/ 2 (1645) — > a 2 7r and aovr and 
their ratio of intensities. There is also reasonable agreement for the line-shape of 7/ 2 (1870) — > 
/ 2 (1270)?7. So masses and widths of these states agree well between the two experiments. The 
evidence for the controversial / 2 (1870) — > a 2 7r signal rests on a small bump in the a 2 7r mass 
spectrum in WA102 data. It now appears likely that some or all of this bump is instead due to 
/ 2 (1910). This state has a strong decay to a 2 7r and weaker decay to / 2 (1270)r?. An earlier WA102 
publication in fact claimed to observe J p = 2 + peaks near 1900 MeV in central production of 
a 2 (1320)7r — > pirn and in / 2 (1270)7T7t [12]. Section 6 presents evidence that t/ 2 (1870) has resonant 
phase variation. Section 7 summarises results and draws conclusions. 

2 Methodology and formulae for fitting data 

The data are fitted by the maximum likelihood method, i.e. fitting every individual event 
without binning. Log likelihood is normalised so that a change of 0.5 corresponds to a change 
in x 2 of 1. For the high statistics available here, log likelihood follows the x 2 distribution closely 
as the number of variables is varied. The following channels are fitted: 
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Here a stands for the inr S-wave amplitude. An incoherent phase space background is also 
included. It arises from experimental cross-talk between the 7/3tt final state and other final 
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states, e.g. Att° and rjr)7r°7r°. This background is known accurately and is discussed in the earlier 
publication [3]; it is in the range 7.7-9.6%, increasing slowly with beam momentum. 

We can dispose of channels (4) and (5) quickly. Their tuti] peaks are narrow and have no 
significant impact on other channels. The /i(1285) is fitted as decaying purely to a^ir. Its mass 
and width need tuning by a few MeV to fit the height and width of the observed peak. The 
77(1440) is fitted with decays to a^ir and 77a, with interference between them. There is also 
evidence for 77(1440) — > /o(980)t/, discussed in Ref. [13] and confirmed in [T4]. The 7/(1440) 
appears clearly only at low beam momenta up to 1200 MeV/c and the fit to it does not change 
significantly from that reported earlier. 

Tests have been made for additional resonances produced in pp — > X + 71 where X has 
quantum numbers J p = 0~, 1 + , 3 + or 4 + . There is no significant evidence for any of these. 
The high spin states would be easily detectable from their strong angular dependence. Adding 
0~ — >■ a vr and f (980)7/ does give a small improvement in log likelihood, typically 20, but this is 
because these amplitudes have no angular dependence and are prone to picking up noise; there 
is no discernable optimum in log likelihood as 777777 mass is varied. The final fit omits 0~ states 
other than 7/(1440). 

2.1 Treatment of partial waves 

Most channels involve a two-stage process pp — > X + tt, X — > Y + n. The orbital angular 
momentum in the decay to X + 77 will be denoted by I and that in the subsequent decay to 
Y + 77 by L. The pp initial state is a mixture of spin singlet and triplet partial waves. The total 
spin S is limited to S z = ±1 or along the beam direction; S z = ±1 give identical angular 
distributions. A problem in the analysis is the absence of polarisation data. The consequence 
is that triplet contributions to X are not cleanly separated between J = £, £ ± 1 and i ± 2. It 
is therefore not possible to do a full partial wave analysis of both production and decay. This 
would become possible if data were available at a future date from a polarised target. 



Fig. 2 sketches the process of production and decay. The 02 and recoil pion are drawn in the 
pp rest frame. The decay of the a<i is shown after a Lorentz transformation to the ci2 rest frame. 
One way of writing amplitudes is to use rotation matrices to express the initial state \J,J Z > in 
terms of a linear combination quantised along the direction. This combination is invariant 
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Figure 2: A sketch of pp — > <22(1320) — > rjix. 
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under a Lorentz boost to the rest frame of the a 2 [T5]- Then the a 2 decay amplitude may be 
expressed in terms of the usual Legendre polynomials. The problem with this approach is that 
the rotation matrices at the first step depend on J. 

This problem may be avoided by a procedure known as the Wick rotation. After a Lorentz 
boost to the rest frame of the a 2 , a rotation of axes is made through an angle —9, the same as 
in the production process but with opposite sign. Rotation matrices then cancel between the 
first and second steps. The Wick rotation alters the angles of tc and 77 in the pp rest frame, 
because of the Lorentz boost. It preserves the fact that the initial state is restricted to J z = 
or 1. Since J is not known, it is however necessary to discard the angular dependence of the 
production process or parametrise it empirically. The way the programme is written, the Wick 
rotation is made for every channel and every event just once, and the required amplitudes are 
stored. This speeds up the analysis by a large factor. 

Two improvements of details over the earlier analysis are introduced. For singlet states, 

S z = 0. The processes pp — > X + tt may go via emission of a pion with orbital angular momentum 

£ > 1, because of the pseudoscalar nature of the pion. For low momenta of the spectator pion, 

several channels rise steeply, see Figs. 7 and 8 below. This is consistent with P-state production. 

If the 77 2 7r final states in reactions (7)-(9) are produced via P-state pion emission, the initial 

state is restricted to J PC = 2~ + unless it is exotic (J PC = 3~ + or 1~ + ). It would be surprising 

if exotics couple to pp and there is no evidence for such exotics in other CBAR data in flight. 

This leads to the useful restriction that the initial state is spin singlet, with S z — 0. A further 

point is that pp — > J PC = 2 h — > [2 h + 7t]l=i has Clebsch-Gordan coefficients such that the 

final state is purely \ J' = 2, J' z = ±1,L — 1,L Z — =Fl >, where J' is the spin of the r] 2 . This 

leads to a distinctive angular dependence for the whole amplitude describing both production 

and decay. It is helpful in isolating the process pp — > 7Ti(2245) —> [t72(2030)7t]l = i. However, 

in addition we detect some significant production from initial spin triplet states, particularly 
PpJ pc = 2++ ^ [2 - + + 7t]l=o 

Interferences between all channels are included. However, several spin triplet states produced 
from pp may feed a single final state such as [J2^]l=i- As a result, interferences between channels 
are not fully coherent. To accomodate this detail, each interference term is multiplied by a 
coherence factor which is allowed to optimise in the range ±2. There are also interferences 
between two a 2 7r and two a 7r combinations for each resonance in rjinr. These interferences are 
fully coherent for a single resonance. 



2.2 The treatment of phase space 

The 772(1870) lies close the the / 2 (1270)?7 threshold. The intensity of the fir\ decay needs to 
be parametrised so as to include the line-shape of the / 2 into the available phase space. The 
formula for the general case pp — > X + Z, where X and Z both have significant width, is given 
by Eq. (40) in Ref. [TB]. This formula is used for channels (2) and (3), pp — > / 2 (1270)a and 
a 2 cr. For the simpler case of ^ 2 (1870) — > / 2 (1270)?7, it reduces to 

(f \ f iV ~ s - m ^ ds! A\p\ MT( Sl ) 

J Ami 7T ^/SS{{M Z — S\) Z + {MT{Si)Y 



where p is the momentum of the 77 in the / 2 ^ rest frame; -y/s and y/s\ are the corresponding 
masses of / 2 ?7 and / 2 . Also FF(s) is a form factor for ?7 2 (1870) — > / 2 (1270)r]. It is taken as a 
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Gaussian exp(— ap 2 ), where a = 4.5 (GeV/c)~ 2 , corresponding to a radius of interaction 0.73 
fm for the overlap of / 2 and rj. From a wide range of CBAR and other data, a is known with 
an error of ±1.0 (GeV/c)~ 2 . In the range of the present data, results vary little over the range 
a = 3.5 — 5.5 (GeV/c)~ 2 . However, the exponential dependence may be an approximation. 




1.5 1.6 1.7 1.8 1.9 2 2.1 

Mass(GeV) 

Figure 3: /2 (1270)7/ phase space (full curve) and an approximation with a Fermi function 
(dashed). 



Fig. 3 shows / 2 77 phase space v. mass. It peaks at 1.96 ± 0.03 GeV, and rises through half- 
height at 1.80 GeV. It is desirable to include this s-dependence into the line-shape of t/ 2 (1870) at 
least approximately. This is done by approximating the phase space of f^t] by a Fermi function 
shown by the dashed curve of Fig. 3: 

F(s) oc 1.0/(1.0 + 1.4exp(5.11(1.76 2 - s))). (14) 

Then the Breit-Wigner amplitude for 7/2(1870) — > / 2 ?/ is 

J M 2 -s-iM[T l + T 2 F{s)Y K ' 

where Ti and T 2 are constants describing decays to (1) a 2 7r and a 7r, (2) f 2 rj. Above 1.96 GeV, 
the Fermi function may be an approximation, but the line-shape of t/ 2 (1870) is falling fast there. 
In principle, this could lead to ambiguities in fitting the t/ 2 (2030), but we find that this state is 
produced in a different range of beam momenta, so in practice there is no problem. 

A further potential complication is that (M 2 — s) of the Breit-Wigner denominator should 
strictly be replaced by (M 2 — s — m(s)) with 

[S) vr J (s> - s)(s> - M 2 )' 1 j 

m(s) is the 'running mass, which makes the formula fully analytic [17J. At a sharp threshold, 
m(s) peaks strongly at the threshold. However, for a threshold as wide as / 2 ?/, its effect is small 
and can be absorbed into optimised values of M and T 2 . 

In our earlier publication on the present data in the year 2000, a sizable ao(980)a amplitude 
was fitted. It accounted for ~ 15% of the 7/371"° cross section, see Fig. 2(a) of Ref. pj. It was 
fitted using the 1993 parametrisation of the a amplitude by Zou and Bugg [H]. Since then, the 
understanding of the a amplitude has improved greatly through the work of Caprini et al. [19J, 
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Figure 4: Mass spectra for (a) 3tt, (b) rjim, (c) tjtt and (d) irn at a beam momentum of 1642 
MeV/c. Histograms show fits to data. 



using the Roy equations. That work fits the tttt scattering length to (0.220±0.005)m~ 1 and alters 
the s-dependence of the tttt amplitude below 600 MeV quite significantly. In addition, there is 
now evidence for significant coupling of the a to KK: g 2 KK (a) / g^(a) = 0.6 ±0.15 [20]. We now 
use Eqs. (l)-(ll) of [2TJ for the ex; they fit accurately the predictions of Colangelo et al. up to 1 
GeV and fit better the inelasticity required above the KK threshold. The additional structure 
in the tttt amplitude is modest, but enough that the evidence for the ao(980)a amplitude in 
present data almost disappears. The remaining signal is barely a 2 standard deviation effect. It 
is now omitted and systematic errors covering the possible signal will be included in branching 
ratios discussed below in subsection 4.2. The earlier ao(980)cr signal was perturbing the fit via 
interferences with ao(980) signals from decays of 7/2(1645), 7/2(1870) and r/ 2 (2030). A consequence 
is that there are now rather large changes to the branching ratios of the 7/2(1870) and r/2(2030) 
to ao(980)7r. The basic difficulty here is that the earlier broad a amplitude had little structure 
and gave a rather flexible fit to the data. The branching ratios of 7/2(1645) and 77 2 (1870) to 
/2(1270)r/ and ci2(1320)7r are more robust and change little. 

The a — > amplitude is needed only up to ~ 1100 MeV. Its form for elastic scattering is 
known quite precisely. In elastic scattering, it is parametrised in the form N(s)/ D(s), where 
the numerator contains an Adler zero just below threshold, making the amplitude weak at low 
momenta. In some production processes with large momentum transfers, e.g. J/^f — > ua, the 
numerator needs to be replaced with a constant in order to reproduce a broad peak in the mass 
range 450-500 MeV, produced by a pole in D(s). In the present data, there is no evidence for 
this behaviour, so the amplitude is taken to be that of elastic scattering. 
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3 Fits to data 



In order to demonstrate the existence of 772 (1645), 772 (1870) and /7 2 (2030), it is necessary to fit 
data at all beam momenta and show that various selections of events require the presence of all 
three, with consistent masses, widths and ratios of decay amplitudes at all beam momenta. The 
earliest publication studied data at just two beam momenta, 1200 and 1940 MeV/c [2J. It was 
immediately obvious that these two momenta required two resonances 772(1645) and 7/2(1870), 
produced with considerably different relative intensities at the two momenta. The 7/2(1870) has 
a strong decay to / 2 (1270)?7 and the 77 2 (1645) does not. The picture developed further when data 
at all nine beam momenta were available. A further 7/2 (2030) was required, with a distinctive 
decay to [a 2 7i}L=2 Here we shall not repeat this lengthy story, but refer the reader to the 
original publications. The conclusion from all beam momenta combined is that just these three 
states are sufficient to fit all the data, with consistent decay amplitudes at all beam momenta. 
Here we simply illustrate the quality of fits to mass projections. 

Fig. 4 shows mass projections for 3tc, mcr], tttj and tttc at a beam momentum of 1642 MeV/c. 
Points with errors are data; fits are shown as histograms. In (a), there is a high mass peak due 
to 7r 2 (1670)77, 7r 2 — > / 2 (1270)7r. In remaining panels there are peaks due to /i(1285), ao(980) 
and a 2 (1320) and a shoulder due to /2(1270). 

At this beam momentum, the irirr] mass spectrum does not distinguish 772 (1645), 772 (1870) 
and 77 2 (2030) cleanly. It is necessary to select events in the mass range 1500-1750 MeV to study 
properties of 77 2 (1645). Figs. 5(a) and (b) show tti] and 7nr mass spectra for this selection. Panels 
(c) and (d) show mass projections for the imrj mass range 1775-1975 MeV, centred on ?7 2 (1870); 
Figs. 6(a) and (b) show projections for the tiht] mass range 2.0-2.25 GeV. Taken together 
with angular dependence in the data, these projections constrain fits to 77 2 (1645), 7/ 2 (1870) and 
7/2(2030) and their individual decay modes. 
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750 1000 1250 1500 250 500 750 1000 1250 
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1000 1500 2000 500 1000 1500 
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Figure 5: Mass spectra at 1642 MeV/c for (a) and (c) 7r°77 in two ranges of M(7nrr)); (b) and 
(d) 7r°7r°. Histograms show fits to data. 
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Figure 6: Mass spectra at 1642 MeV/c for (a) and (b) 7r°r/ and 7r°7r° in the M(7r7r7/) mass range 
2.0-2.25 GeV; (c) the nnr] mass spectrum selecting a 2 and vetoing / 2 (1270) (d) the converse 
selection. Histograms show fits to data. 



In order to display more clearly the 772(1870) , Figs. 6(c) and (d) shows mass projections 
which select a 2 (1320) and reject / 2 (1270) or vice versa. In (c), a 2 7r events are selected with 
M(irri) in the mass range 1318 ±55 MeV; events with 7r°7r° in the mass range 1.0-1.455 GeV are 
rejected. In (d), / 2 (1270) is selected in the mass range 1.1-1.37 GeV and events containing 7rr) 
in the mass range 1.15-1.435 GeV are vetoed. The strong / 2 (1870) peak in / 2 (1270)?7 is clearly 
visible in Fig. 6(d). In Fig. 6(c), one sees the combined rj7nr mass spectrum from r/ 2 (1645) and 
772(1870) — > a 2 7r. Similar cuts in data selection illustrate the presence of ?? 2 (2030) — > [a 2 7r]i =2 . 
Fits to data are of similar quality at all beam momenta; further examples were shown in Ref. 

®- 

In order to achieve good fits to data, the eleven channels of Section 1 need to be sub-divided 
to include separate channels for (i) 77 2 (1645) —> a 2 7r and a n, (ii) 77 2 (1870) — > / 2 (1270)?7, a 2 7r 
and a vr, (iii) ?7 2 (2030) ->■ [a 2 vr] i=2 , [a 2 7r] L=0 , a ir and / 2 (1270)t/, (iv) / 2 (1910) -> [a 2 7r] i= i and 
[/ 2 (1270)77] L=1 , (v) / 2 (2001) -> [a 2 7r] L=3 only, (vi) / 2 (2240) -)• [a 2 7r] L=1 , 3 and [/ 2 (1270)7r] L=1)3 . 
At each beam momentum, phases for all these channels are fitted freely. 

In earlier work, there was some evidence for r/ 2 (1870) — > [i]a]L=2- There remains a small 
improvement in log likelihood when this channel is included. This improvement is however 
typically 30, which is less than for almost all other channels. Furthermore, its inclusion increases 
errors on other channels, i.e. fits become less stable including it. The problem is that there is 
no narrow signature of the o — > tttt S-wave, so it tends to absorb any noise in the data. This 
decay is now omitted. 



4 Results 

Table 1 shows changes in log likelihood when each channel is omitted from the fit and all other 
channels are re-optimised. One sees immediately the significance level of each channel and its 
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Table 1: Changes in log likelihood when each channel is removed from the fit and others are 
re-optimised. 



dependence on beam momentum. 

We now wish to draw conclusions from peaks observed in some channels of data. Fig. 7 
shows cross sections v. beam momentum. The absolute normalisation is taken from cross 
sections determined in Ref. [3] for the whole r]3iT data; they are uncorrected for branching 
fractions of 77 and 7r° — > 77 and therefore correspond directly to the number of events collected. 
The integrated cross section varies little with beam momentum, with a small (10%) enhancement 
near 2270 MeV; it is shown in Fig. 2 of Ref. [3]. On Fig. 7, smooth curves for individual channels 
(2)-(12) are drawn through the data within one standard deviation. Errors are typically 5 — 10% 
and some examples will be displayed with errors in Fig. 8. Cross sections for several channels 
expand at low beam momenta 900 and 600 MeV/c to follow the \jv dependence of the total 
cross section, where v is the beam velocity in the centre of mass. 

In Fig. 7(b), there is a strong peak in the channel f 2 ao near 2250 MeV and a broader peak 
in 772(2030)77 at slightly higher mass. There is a possible peak in 7/2(1645)7? in the same mass 
range. There is also a peak in / 2 (1910)7r near 2150 MeV. The channel 7t 2 (1670)t7 peaks strongly 
at high mass. There is evidence for production of / 2 (2240)tt at the highest two beam momenta. 
Other channels show only weak structure, except that a 2 cr peaks at low masses. 

The question arises how to interpret the peaks. Are they due to resonances? Here it is 
necessary to take care over details in the formulae. The amplitude for a resonance such as 
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Figure 7: Cross sections for each final state. Smooth curves are drawn through cross sections 
for each channel. 



a 2 (2255) / 2 (1270)a is 

A = ^p(pp, s)p(f 2 a ,s)/[k(M 2 - s - iMT tot )]. (17) 

The factor k in the denominator is the momentum of the p in the pp centre of mass; this is 
the flux factor for the incident beam. The phase space factor for p(pp, s) in the numerator is 
k/y/s, multiplied by a centrifugal barrier factor. We choose to make comparisons with data by 
accounting explicitly for the intensity factors kj i/i for pp and 1/k 2 from the flux factor; the two 
together give a factor l/k\fs = 1/v, where v is relativistic velocity. We multiply cross sections of 
Fig. 7 by a factor v/vq, where v is evaluated at 2410 MeV, the highest data point. Remaining 
centrifugal barrier factors and form factors are less certain and are modelled in the fit to results. 

Fig. 8 shows data for (v/vo)a. For /^o m Fig. 8(a), there is a peak closely resembling 
the a2(2255), reported in the earlier combined analysis of 7 = 1, C — +1 CBAR data with 
M = 2255 ± 20 MeV and T = 230 ±15 MeV; it was observed as a clear peak in data for 
pp — > / 2 (1270)7r, see panel (£) of Fig. 38 of Ref. [5]. It was observed in both 3 F 2 and 3 P 2 decays 
to / 2 (1270)7r, with an amplitude ratio —2.13 ±0.20 favouring coupling of pp to 3 F 2 . The dashed 
curve of Fig. 8(a) shows the remaining phase space factor calculated with this ratio for 3 F 2 and 
3 P 2 production and with S-wave decay to / 2 ao- The rise of this curve with mass fails to fit the 
data. Any other ratio of 3 P 2 and 3 P 2 also fails to fit the data. Adding P-state / 2 ao decays makes 
the dashed curve peak even higher in mass. 

The full curve shows the result of multiplying the dashed curve by the line-shape of a 2 (2255), 
assuming a Breit-Wigner resonance of constant width. It is remarkably close to the data. We 
regard this as further confirmation of the a 2 (2255). A small improvement is possible by reducing 
the mass by 10 MeV to 2245 MeV and decreasing the width to 225 MeV, but these changes are 
within the errors of the earlier determination and also within errors of present data. 



12 




2000 2200 2400 2000 2200 2400 

Mass (MeV) 



Figure 8: Production of (a) / 2 (1270)a compared to phase space for S-wave production (dashed 
curve), a fit to j 2(2255) (full curve) and a fit with parameters of this state optimised (dotted 
curve); (b) r/ 2 (2030)7r compared to phase space for P-state production (dashed curve), a fit using 
the line-shape of 7r 2 (2245) (full curve) and a fit with parameters of this state optimised (dotted). 



Fig. 8(b) shows results for the channel ?72(2030)7r. Production of this final state is dominantly 
(~ 74%) with pp helicity 0, suggesting production via pp — > 7r 2 (2245): 

pp ->• tt 2 (2245) ->■ [77 2 (2030)tt] £=1 (18) 
^(2030) [a 2 7t] L=2 (19) 
a 2 (1320) -> [r]7i] L=2 . (20) 

The full amplitude D for this process is rather distinctive. Suppose it is written in fixed axes in 
the pp centre of mass. The initial state is spin singlet, helicity 0. Suppose the spectator pion of 
Eq. (18) is produced at angle 9 of Fig. 2 and with azimuthal angle <fi around the beam direction. 
Let the momentum of the decay pion in Eq. (18) be q. Let the decay angle of the a 2 (1320) in 
Eq. (20) after the Wick rotation be a, with azimuthal angle (3, using the same axes as for Eq. 
18. Then 

D oc J^qBtiq) sin 8 sin a cos aie^-® - e - *^) (21) 

= iy/2qBi(q) sin 6* sin a cos a sin (0 — (5). (22) 

Here Bi(q) is the centrifugal barrier factor for L = 1 decay. The amplitude for the decay 
a 2 (2030) — > / 2 (1270)?7 has an identical form. Amplitudes for decays of a 2 (2030) — > [a (980)7r] i=2 
take a similar form with angles a and those of the decay pion. The amplitude for decays to 
a 2 (2030) [a 2 (1320)7r]i =2 involves a combination of spin 2 of the a 2 (1320) with L = 2 to make 
spin 2 of the a 2 (2030). This amplitude may be constructed along the same lines using Clebsch- 
Gordan coefficients or using tensor algebra. Because of the additional angular dependence on 
L = 2 in the decay to [a 2 (1320)7r]i =2 , it is particulary distinctive and stands out clearly in the 
data. 

From earlier work, the parameters of 7r 2 (2245) have sizable errors, M = 2245 ± 60 MeV, 
T = 320^40° MeV The full curve of Fig. 8(b) shows a fit with these parameters; the width 
is clearly too large and the mass somewhat too low. 



13 



Production of 7/ 2 (2030)7r is (26 ±4)% in intensity from pp helicity 1 at the peak; this helicity 
1 component may be accomodated with a,2(2255) — > S-wave 7/2(2030)71". Note that a2(2255) 
cannot contribute to [r/2(2030)7r]^ = i. Also [7/2(2030)7r]^=2 has a strong variation with mass and 
will be strongly suppressed by the I = 2 centrifugal barrier. The dotted curve of Fig. 8(b) 
is obtained by adjusting the mass and width of 7r 2 to M = 2285 ± 20 (stat) ± 25(syst) MeV, 
T = 250 ± 20(stat) ± 2h{syst) MeV and including 26% of the intensity via a 2 (2255). 

The data determine relative phases of strong channels with errors of ~ ±15° at every mo- 
mentum. The phase variation of /2&0 and a2(2030)7r channels agree within these errors. If f^ao 
is resonant, then so is 02 (2030) tt and this can only be explained in terms of some combination 
of 772(2245) and a2(2255). There are further triplet states as(2275) and ai(2245) in this vicinity, 
but they can only contribute to decays to [i]2(2030)7i]e=2', this is inconsistent with the observed 
line-shape of the peak and can only make a small contribution to the data. Systematic errors 
assigned to the mass and width of 7T2(2245) cover any weak contribution from this source. 

4.1 Other channels 

The production of 7/2(1875)7r is distinctively different to that of 772(2030)77. Production from 
the initial pp system with helicity ±1 is strongly dominant, requiring an initial triplet state. 
However, no conclusion can be drawn from the slowly varying cross section in Fig. 7(a). A minor 
correction to Ref. [3] is that the 7/ 2 (1875)7r signal plotted there in Fig. 2(a) was multiplied by a 
factor 2 to make it clearly visible; the intensity recorded there is close to the present analysis. 

Production of / 2 (1910)7r is roughly equal from initial states with helicity and 1, but shows 
significant slow variation with mass. This is possible from initial triplet states. The peak near 
2150 MeV may indicate production from the initial state a2(2175). 

Production of 7t 2 (1670)t7 is dominantly via pp helicity 1, (i.e. spin triplet) but again fluctuates 
smoothly with mass by more than a factor 2. No firm conclusion can be drawn from the variation 
of cross section with mass. 

A final point is that data are available at 900 MeV/c and were included in the earlier anal- 
ysis [3]. However, there are small cross sections for several channels and these are difficult to 
determine with confidence. This beam momentum has been studied, but is discarded from the 
present analysis because of substantial systematic errors in weak channels. There are also low 
statistics at 600 MeV/c, and this momentum is omitted for the same reason. 

4.2 Branching fractions 

Fortunately, the branching fractions of 77 2 (1645), 7i 2 (1870) and 7/2(2030) are not sensitive to the 
question whether or not their production goes directly via resonances in the pp channel. If such 
resonances are involved, all decay modes of these channels pick up the same phase variation 
from the production process. The dominant decay of r/2(2030) is to [ci27r]L=2; that for 7/ 2 (1870) 
is to /2(1270)?7 and that for 7/2(1645) is to a27r. Branching fractions will be quoted with respect 
to these dominant channels. Some of the smaller branching fractions have changed significantly 
since the analysis of the year 2000. 

A general comment is that data determine ratios of amplitudes rather than ratios of intensites. 
That is, log likelihood has a parabolic minimum as a function of the amplitude ratio. The 
procedure is therefore to fit data using, for example, the magnitude of the 7/2(1645) — > a 2 TT 
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amplitude as one variable and the ratio of amplitudes 77 2 (1645) — > ao 7r / ? ?2(1645) — > a 2 n as 
the second. This makes it easy to investigate fluctations with beam momentum. The average 
of each ratio of amplitudes is found by weighting values at each beam momentum with the 
total number of a 2 (1645) events summed over both channels. The same procedure is used for 
77 2 (1870), 772(2030) and f 2 Ti channels like f 2 (l9l0)ir. This gives more stable results than trying to 
determine errors of ratios at each momentum and using them in the weighting procedure; errors 
fluctuate somewhat from momentum to momentum because of large errors for phases fitted to 
weak channels. Having found the average ratios of amplitudes, a second pass is made through 
the fit, fixing these ratios at all momenta. Then the ratio of intensities is obtained from the 
total number of events fitted to each channel, integrated over all beam momenta. 

In evaluating branching ratios, it is necessary to make use of Eqs. (13)-(15). Here there is 
a dilemma. It is convenient to parametrise the Breit-Wigner denominator as far as possible 
with constant widths for each decay channel, as is conventional in the Particle Data Tables. 
Our procedure is to modify the numerator of Eq. (15) to T 2 F(s), so as to agree with the 
denominator. We then evaluate branching ratios at each beam momentum, and average over 
momenta to determine Ti/T 2 . This procedure converges within errors after one iteration. The 
same procedure is used to evaluate the effects of L = 2 centrifugal barriers on widths for 
r7 2 (2030) ->■ [a 2 ir] L = 2 and 77 2 (2030) ->• [a 7r] L=2 . 

The acceptance for strong channels / 2 (1270)?7 and [a 2 7r]i =2 varies rather strongly with beam 
momentum. Numerical results are shown in Table 2, normalised to 1 at the highest beam 
momentum 1940 MeV/c; they are evaluated using the Monte Carlo simulation of the detector. 
One sees a large variation of the intensity ratio R for some channels. This variation almost 
disappears at 1940 MeV/c, where all decay channels are nearly fully open. We shall tabulate 
branching ratios R corrected to this momentum. Above this momentum, results may be affected 
by errors in the Fermi function adopted in Eq. (14), so this is close to the optimum compromise. 
Errors in branching ratios due to these uncertainties are included in errors quoted in Table 3 
below. We regard this procedure as an improvment on the work of Ref. [3], where branching 
ratios were evaluated purely from geometric acceptance without the effects of centrifugal barriers 
or the form factor FF(s) of Eq. (13). 



moment um( Me V /c) : 
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a ir/ f 2 rj 


2.01 
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1.0 


V2 (2030) 


[a 2 7r] L=0 /[a 2 7r] L=2 


1.78 
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1.0 




a n/[a 2 7v] L=2 


1.57 


1.43 


1.26 


1.15 


1.09 


1.04 


1.0 




f 2 r]/[a 2 n] L=2 


0.98 


1.00 


1.01 


1.02 


1.01 


1.01 


1.0 


/ 2 (1910) 


f2V/a 2 n 


0.738 


0.803 


0.871 


0.917 


0.948 


0.976 


1.0 



Table 2: Variation of acceptance (including form factors) with beam momentum 



Coming to technicalities, it will be necessary to correct the number of observed a 2 (1320)7r 
and / 2 (1270)?7 events in 7r°7r°r/ for unobserved decays. The / 2 (1270) has a branching ratio of 



15 



0. 848/3 to 7r°7r° [TU]. The amplitude for the 1 = 1 component of pp — > ttX, X — > f 2 rj is given 

by 

A(I = 1) = 7r?X 2 ° 3 + + rr° 3 X° 32 (23) 

In (23), X stands for the amplitude of X — > f 2 f]; the brackets in the last line identify pions 
coming from / 2 (1270). For X — > a 2 n, the decay amplitude of X is y / T/3(7r + a 2 — n ^ + n~a 2 ), 
where the minus sign for ir°a 2 can be absorbed into the fitted phase for this channel. In the 
TT°a 2 final state, what is actually observed is 

b(i = i) = + (4vK\ + AMvX + + A(4v)4 + , 

(25) 

1. e. six combinations. The coherent sum of all combinations is fitted to the data. 

Table 3 shows in column 3 branching fractions from the previous analysis [3] for comparison 
purposes. The next column lists what is fitted now. The final column corrects this for all charge 
states and for the branching fractions of 02(1320) — > nr] (14.5%) and a (980)7r to nrj. This last 
branching fraction is taken to be the value used by WA102, (86%) for easy comparison with their 
results. This is close to the value adopted by the PDG [10J. If their value is adopted, ao(980) 
branching fractions increase by a factor 1.015 ± 0.021. Results in columns 4 and 5 supercede 
the earlier results. 

Branching ratios of decays to ao(980)7r final states are unstable. They depend somewhat on 
whether or not the a (980)cr channel is included in the fit. The basic difficulty is that the a 
amplitude is broad and gives rise to interferences all over the 4-body phase space. The a Q a 
channel is therefore not well determined. In the absence of definite evidence that it is needed, 
we omit it. All signals are weak and their phases with respect to dominant decays have quite 
large errors. There are also strong correlations between couplings of 7/2(1645), 7/2(1870) and 
?72(2030) to aovr. In the earlier analysis, aovr decays interfered with ao<j, giving apparently small 
but unreliable errors. 

From the present analysis, the ratio for 77 2 (1645) lies close to the WA102 result 0.077 ±0.016; 
the ao(980) signal is clearly visible as a peak in their raw data [4]. For this reason, the weighted 
mean of these two values has been adopted and fixed for r7 2 (1645) . This helps stabilise fits to 
a n. However, there are large correlations between 77 2 (1870) and ^ 2 (2030) decays to a (980)7r. 
Together with uncertainties due to possible contributions from ao(980)a, the result is that 
branching fractions of 7/2(1870) — > and ?72(2030)7r can both vary freely over the range 0.1 
to 0.85 in present data. Errors on values given in Ref. [3] for these two channels need to be 
increased substantially to take account of these systematic errors. New estimates are given 
in Table 3 . These uncertainties are not significantly correlated with branching fractions to 
/ 2 (1270)?7 and a 2 7i. Branching ratios between / 2 (1270)?7 and a 2 (1320)7r decay modes are mostly 
quite stable, because these decays are strong. The ratio of decays of 77 2 (1870) to / 2 (1270)r/ and 
a2(1320)7r remains stable. So does the ratio of decays of 7/2(2030) to f 2 (1270)r] and [a 2 T[}i=2, 
because these are conspicuous decays. 

A mistake has been located in the branching fraction of r/ 2 (2030) to [a27r]x, = o reported in Ref. 
[3]. The corrected value is given in entry 6 of Table 3. This decay is much weaker than that to 
[a27r]L =2 . 
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0.43 ±0.15 


0.75 ±0.34 


0.13 ±0.06 


772(2030) 


Si2[a 2 7r] L= o/Bi2([a 2 7r] L=2 ) 




0.05 ±0.03 


0.05 ±0.03 


/ 2 (1910) 


BR{f 2 r,)/BR{a 2 Ti) 




0.54 ±0.29 


0.09 ±0.05 



Table 3: Branching ratios; column 4 shows values for present data and column 5 shows values 
corrected for all charges and all decay modes of / 2 , a 2 and do- 



Beam momentum (MeV/c) 


amplitude ratio r/ 2 (1870) — > a 2 7r/r/2(1870) — > / 2 r7 


1050 


0.383 ±0.067 


1200 


0.321 ±0.034 


1350 


0.219 ±0.050 


1525 


0.202 ±0.058 


1642 


0.320 ±0.050 


1800 


0.334 ±0.063 


1940 


0.308 ±0.065 



Table 4: The ratio of amplitudes ^(1870) -> a 2 (1320)7r/r/ 2 (1870) ->■ / 2 (1270)r/ at individual 
beam momenta. 



Table 4 gives the important branching ratio of amplitudes r/ 2 (1870) — > a 2 (1320)7r/r/2(1870) — > 
/ 2 (1270)r/ at all beam momenta, in order to illustrate the stability. We choose to take the 
simple mean 0.298 over beam momenta, so as to avoid bias from fortuitous fluctations in errors 
with beam momentum. Fluctuations about the mean are above statistics by a factor 1.4. The 
statistical error is increased to allow for this in Table 3. 

In earlier work, the error on the mean was taken as the statistical error divided by y/~N, 
where N is the number of beam momenta. However, it is now clear that systematic errors are 
somewhat larger than this. The systematic errors are estimated from (i) variation of results 
with the ingredients included in the fit, particularly the number of interferences included be- 
tween channels; (ii) variations with masses and widths of 77 2 (1645), 772 (1870) , r/ 2 (2030), / 2 (1920), 
/ 2 (2001) and / 2 (2240), (iii) possible contribution from / 2 (2293), (iv) uncertainties in the back- 
ground from other final states. Ultimately, systematic errors dominate for all branching ratios, 
particularly for decays to the weak a^ir channels. 

5 Fits to WA102 data 

The WA102 collaboration measured central production of ttttt] and produced separate sets of 
data for 77 — > 77 and 7r + 7r~7r° |4]. These data have been read from their graphs and refitted. 
Their approach was to use the K-matrix for t/ 2 (1645) and t/ 2 (1870). If a single amplitude is 
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mass ( GeV) mass ( GeV) 



Figure 9: Fits to WA102 data for central production of a^i^-, cto 71 " an d /2(1270)?y with J PC = 2 h . 
The first column shows data for 77 — > 77 and the second data for 77 — > 7t + 7t~tt°. Full curves are 
fits without ?72(2030) and dashed curves with it. 



used in this approach, the amplitudes for 77 2 (1645) and 77 2 (1870) each loop round the Argand 
circle once; continuity of the amplitude then requires a zero between them. The dip due to this 
zero played a significant part in fitting the data. 

The basic assumption of the K-matrix approach is that resonances combine in a production 
process in an identical way to elastic scattering. This is a pure assumption. In elastic scattering, 
ingoing and outgoing waves for all coupled channels must sum asymptotically to unit intensity. 
However, in production processes considered here, the amplitude is only a tiny fraction of the pp 
total cross section. There is no obvious reason why 2-body unitarity should apply in the same 
way as for elastic scattering. The K-matrix approach has been tested on four sets of experimental 
data in Ref. [22]; it failed seriously in every case. If it is used, phases of r? 2 (1645) and ?72(1870) 
need to be fitted freely, since final states may rescatter between one another, generating phases 
which need to be fitted arbitrarily. Our analysis therefore uses the isobar model. This makes no 
assumption about the effects of unitarity and allows separate phases for 77 2 (1645) and 77 2 (1870). 

Refitted results are displayed on Fig. 9. Data for 77 — >• 77 and 77 — > 7r + 7r~7r have been fitted 
simultaneously, including decays to all of ao(980)7r, / 2 (1870)r7 and a2(1320)7r. They have been 
fitted with and without 772 (2030). That component was not known at the time of the WA102 
analysis. Full curves on Fig. 9 show fits without 7/2 (2030) and dashed curves the fits including 
it. With it, the total \ 2 improves from 124.5 to 95.8. The reason is obvious: it provides extra 
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freedom in fitting small defects in the mass region above 2 GeV. This improvement in \ 2 of 28.7 
needs to be balanced against the fact that there are three extra complex coupling constants for 
the three channels of 7/2 (2030), i.e. 6 extra fitting parameters. The improvement in \ 2 is 3.6 
standard deviations. It is debatable whether or not 772 (2030) is really present. Conclusions will 
be drawn here from fits without it; these fits are more secure. 

A significant point is that the K-matrix zero between the 1645 and 1870 MeV peaks made 
both of them narrower and pushed the peaks apart. In the WA102 fit, the mass of 7/2(1645) was 
1605 ± 12 MeV for 77 — > 77 data and 1619 ± 11 MeV for 77 — > 3n. These are to be compared with 
the CBAR determination of 1645±6(stat) ±20(syst) MeV. In the isobar model, the dip between 
the two resonances can be filled in by interferences between them. Table 5 shows fitted masses 
and widths for 7] 2 (1645) and 7/2(1870) for two cases. In the first, column 2, the parameters of 
7^2(1870) are fitted freely. 

The width of the 772 (1870) tends to run away to a large value. In Fig. 9(d), the width of the 
7/2(1870) peak in decays to / 2 (1270)7/ is sensitive to scatter in the points near the peak. The 
second fit (column 3) is made adding to x 2 a contribution given by the CBAR masses and widths 
with statistical and systematic errors combined in quadrature. This extra constraint stabilises 
the fit and gives mass and width for 772 (1870) closer to the CBAR values. In our opinion, the 
third column is the more reliable, bearing in mind that the addition of the 7/2(2030) increases 
the uncertainties from WA102 data even further. 




1.5 1.6 1.7 1.8 1.9 2 2.1 



mass (GeV) 

Figure 10: The line-shape of 7/2(1870) from the combined fit to CBAR and WA102 data. 



The explicit re-arrangement of Eq. (15) for 7/ 2 (1870) ->• f 2 (l270)i] is 

= [0.685F( S )Mr] 1 / 2 
J [M 2 -s-i(0.672 + 0.685F(s))Mr]' 1 ' 

where M and V refer to values for 772(1870) in the last two lines of Table 6; F(s) is given by Eq. 
(14). The CBAR data are fitted using Eqs. (14) and (26). Numerical values in the denominator 
are such that \ f\ 2 integrated over s reproduces the branching fractions for a27i, aovr and ^2(1270)77 
in Table 3. 

The line-shape of 77 2 (1870) is shown in Fig. 10 and is very close to that of a Breit-Wigner 
resonance of constant width. It peaks at 1792 MeV. This is related to the opening of the 
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Unconstrained fit With CBAR constraint CBAR values 



X 


1 1 p. fi 
110. o 


1Z4.0 




M(1645) 


1635 ± 12 


1630 ± 9 


1645 ±6(stai) ± 20{syst) 


r(1645) 


252 ± 36 


225 ± 16 


200 ± 5(stat) ± 25(syst) 


M(1870) 


1833 ± 20 


1829 ± 12 


1825 ± 5(stat) ± 15(syst) 


r(1870) 


332 ± 45 


293 ± 24 


221 ± 20(stat)±§(syst) 



Table 5: Masses and widths fitted to WA102 data without any constraint (column 2) and with 
a penalty function given by errors on CBAR masses and widths (column 3); Column 4 shows 
CBAR values of masses and widths. 



/ 2 (1270)t7 threshold. The pole position is at 1798 ± 20 - z(130 ± 12) MeV. If CBAR values are 
used instead, the imaginary part of the pole position decreases to 109 MeV. 

The remaining question is whether there is really any significant a 2 (1320) signal in WA102 
data due to r/ 2 (1870) -» a 2 (1320)7r. In a 1997 publication [12], WA102 claimed to observe 
772(1645) — > a 2 (1320)7r in central production of 4tt, with a small shoulder at high mass which 
could be ?72(1870) . They reported a strong 2 + signal in a 2 (1320)7r — > pirn at 1900 MeV and a 
broad 2 + peak in / 2 7T7t near 2000 MeV. The integrated 2 + signal was considerably stronger than 
that for J p = 2 h . However, further data reported in the year 2000 on central production of 
47T were interpreted in terms of 2 h — > a 2 (1320)7r — > Att, produced only with J z = ±1 [23]. The 
2 ++ — > / 2 (1270)7T7T was found again but no 2 ++ — > a 2 (1320)7r. They make no comment on why 
this change from the 1997 work occurs. In central production via Pomeron exchange, there is 
no obvious reason why 2 ++ should not be produced with J z = 0, ±1 and ±2. Our view is that 
the data really need to be fitted with all allowed values of J z for both J p = 2 + and 2~. It would 
be valuable if the Compass collaboration could check this point. 

In the CBAR data analysed here, there is clear evidence for / 2 (1910) and / 2 (2001) decaying 
to a 2 (1320)7r and / 2 (1270)r7. The a 2 7r decay dominates. This is readily understood from the 
fact that the L = 1 centrifugal barrier inhibits decay to / 2 (1270)r/. We suggest that the small 
bump claimed by WA102 in a 2 (1320)7r at 1860 MeV is due to / 2 (1910). We have fitted WA102 
data using PDG parameters M = 1903, T = 196 MeV for / 2 (1910) instead of r? 2 (1870). The fit, 
shown in Fig. 11, gives a slightly improved description of the data (by 11 in x 2 ), but we are 
unable to go back to the original data and check the J p analysis. 

H 400 

I 200 




Figure 11: Fits to WA102 data with 7? 2 (1870) -> a 2 7r replaced by / 2 (1910) -> a 2 n: (a) r] 77, 
(b) rj — > 7r + 7r~7T°. 
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It remains an interesting question why 77 2 (1870) has a fairly large / 2 (1270)7/ S-wave decay 
The dispersive term m(s) of Eq. (17) leads to attraction in this channel near threshold [17]. It 
necessarily favours this decay mode. There is an isospin partner 7^(1880) for 7/2(1870). It too 
has a strong decay mode to a 2 (1320)?/ [23] . 

6 Evidence that the 772 (1870) is resonant 

Several checks have been made that the 772 (1870) has resonant phase variation. The first check is 
to remove the phase variation by replacing the amplitude by its modulus and refitting the data. 
At all beam momenta, this leads to a highly significant worsening of log likelihood. Column 2 
of Table 5 shows the changes against beam momentum. The definition of log likelihood is such 
that a change of 0.5 in log likelihood should correspond to a one standard deviation change. 
We have already remarked that fluctuations in Table 4 are a factor 1.4 above statistics. A more 
extensive examination of fluctuations in branching ratios shows they are in some cases up to a 
factor 2 above statistics, and that has already been taken into account in errors quoted in Table 
3. Adopting this as a general rule leads to the conclusion that changes listed in Table 3 can be 
equated to changes in \ 2 - They average to 11.6 standard deviations per momentum, i.e. 28<r in 
total. 



Momentum (MeV/c) 


(i) 


(ii) 


1050 


43 


31 


1200 


81 


98 


1350 


106 


120 


1525 


91 


47 


1642 


95 


55 


1800 


195 


127 


1940 


199 


86 



Table 6: Changes in log likelihood when the amplitudes for 7/2(1870) (i) are replaced by their 
modulus (no phase variation), (ii) use a denominator A — m{s) — iMT(s), see text. 

The second check is that including the dispersive term of Eq. (16) into the Breit-Wigner 
denominator has little effect on log likelihood after minor alterations to fitted mass and width, 
well within errors quoted in Table 6. The dispersive term peaks at the /2(1270)7/ threshold, but 
with a large full width of ~ 300 MeV. Numerically, it is easily absorbed into small shifts of fitted 
parameters. 

There remains the possibility that the associated cusp in the real part of the amplitude could 
explain the phase variation without a resonance. The third check is to replace the Breit-Wigner 
denominator of Eq. (15) by 

D = A-m(s) -iM[Y 1 + T 2 F(s)}. (27) 

This removes the resonance by changing M 2 — s to a constant A; it is the term in s which drives 
the real part of the amplitude through zero on resonance. Fig. 12 shows the Argand diagram 
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of the amplitude fitted to WA102 data. [The vertical scale is no longer limited to 1 at the peak, 
because the imaginary part of the amplitude no longer needs to reach 1 in the absence of a 
resonance.] ^ 

0.6 
0.4 
0.2 



" 0.2 0.4 0.6 0.8 1 

Ref 

Figure 12: The Argand diagram for a fit to WA102 data using Eq. (26). 

The imaginary part of the amplitude is forced by the data to peak at ~ 1810 MeV, as before. 
The real part is positive everywhere. It turns out that the data do not force the constant A 
to go negative; instead it optimises close to +0.5. Then fits to 77377° data are worse than for a 
resonance by the amounts shown in column 3 of Table 6. At the lowest three beam momenta, 
where the 772(1870) is strongest, results are similar to column 2. At higher momenta, the changes 
drop by roughly a factor 2. The reason for this drop has been traced to the fact that the phase 
variation in Fig. 12 adds a degree of freedom in the fit compared with column 2, where there is 
no phase variation at all. However, the significance levels in column 3 still average 9 standard 
devations per momentum. If any phenomenologist wishes to develop a more complete dynamical 
model, the data are publicly available from the authors. Meanwhile, the evidence for resonant 
behaviour appears to be strong. 

7 Conclusions 

This work confirms that the 772(1870) has a branching fraction to /2(1270)77 comparable with 
that to a2(1320)7r, in agreement with our earlier analysis. It is not possible to fit these data with 
the large branching fraction found by WA102. Results of the two experiments agree well for 
the mass and width of 77 2 (1645) and the branching fraction of its decays to a (980)77. They also 
agree quite well for the mass and width of 77 2 (1870) from decays to /2(1270)r/. Small branching 
fractions reported in Table 2, particularly for decays to a 77, have changed significantly from 
earlier values for a complex of reasons which are understood. 

The WA102 collaboration has found evidence for a weak decay mode of 772(1645) — > KKir 
[25] (7% of that for decays to 0277); however, they find no evidence for 772(1870) in the same data. 
It therefore appears unlikely that the 772(1870) is the ss partner of 7/2(1645) and 772(1670). Fig. 
13 shows trajectories of J = 0, C — +1 for several sets of quantum numbers. In (b) 772(1645), 
7/2(2030) and 772(2267) are consistent with the 1 D2 trajectory with the same slope as the others. 

Hybrids with J PC = 2 h are predicted around 1900 MeV by Isgur and Paton [26J and 
Godfrey and Isgur [27j. The interpretation given by many authors, including ourselves [3], is 
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a) n (b) n 

Figure 13: Trajectories of resonances for several quantum numbers of 7 = 0, C = +1 states; 
masses are shown in MeV. 

that 772(1870) and 7T 2 (1880) make a hybrid pair, somewhat higher in mass than the 7r 1 (1600) 
of the PDG, whose average mass is now 1662 MeV. It would be valuable to search for ssg 
partners in J/ty decays at BES 3. Page [28] predicts that hybrids will decay dominantly to 3 P 
qq plus a pion, in agreement with the strongest observed decay modes a 2 7r and / 2 ?7 of 772(1870). 
Predictions for their branching ratio are subject to systematic errors from (i) the effect of the 
dispersive attraction from the / 2 77 threshold and (ii) possible mixing of the hybrid with qq ID 
and 2S states. The phase space alone for a 2 7r is 1.184 ±0.023 times that for f 2 r) after integrating 
over the line-shape of 7/2(1870) and taking account of the form factor exp[— (4.5 ± 1.0)g 2 ]. 

Li and Wang propose that the 7/ 2 (1870) is the n = 2 qq state and 7/(2030) is the n = 3 state 
[29] . However, this would required a trajectory with twice the slope of other J p . 

Afonin has presented an interesting scheme to accomodate known light mesons. Its 
general features are appealing. However, we question the way J PC = 2~ + states are included. 
He includes 77 2 (1645), r/ 2 (1870) and r/ 2 (2250) for 7 = and 7r 2 (1670), 7r 2 (2100) and 7r 2 (2245) 
for 7 = 1. The large mass difference between 77 2 (1870) and 7T2(2100) is unexplained; the mass 
splitting between 7/ 2 (1870) and / 2 (1934) is reversed for 7r 2 (2100) and a 2 (2030). 



Momentum (MeV/c) 


Change in log likelihood 


1200 


408 


1350 


387 


1525 


420 


1642 


536 


1800 


484 


1940 


215 



Table 7: Changes in log likelihood when 77 2 (2030) is removed from the fit to 77371"° data. 

His model conflicts with present data, which require the presence of 77 2 (2030); Table 7 shows 
changes in log likelihood if it is omitted and all other parameters are re-optimised. They are on 
average 1.18 x values in Table 7 for 772 (2030) — > [a 2 7r]L = 2, but they are not as large as the sum 
for all decays of 772 (2030). This is because of correlations between decay channels. Nonetheless 
they are still highly significant. Furthermore, there is independent evidence for 7/2(2030) in 
pp -> f]7i7c [6J. In those data, there is a strong peak with M = 2040 ± 40 MeV, T = 190 ± 40 
MeV in [/ 2 ( 1270)7/] L=0 and a smaller, but still significant peak in [/ 2 (1270)77]i /=2 ; unfortunately 
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the resonant phase variation cannot be checked in those data because there is no other strong 
feature in other singlet states near this mass. 

The vr 2 (1880) is listed by the PDG in four sets of data: (i) in r/r/vr with M = 1880 ± 20 MeV 
by Anisovich et al. [31], (ii) in ir'p ->- r]n + iT-ii-p with M = 2003 ± 88 ± 148 MeV by E852 [32], 
(iii) in n'p uir-Tr°p with M = 1876 ± 11 ± 67 MeV in further E852 data of Lu et al. [33], 
and (iv) in n'p T]T]7r-p with M = 1929 ± 24 ± 18 MeV by Eugenio et al. (E852) [21]. Of 
these, the second one could be 7r 2 (2005). 

A further result from the present analysis is that there is evidence for the presence of 
a2(2255) — > /2(1270)a (980) with parameters close to those of Ref. [TT] . It is the third set of data 
in which it has been observed, the others being irr] and 37r°. There is also evidence that the chan- 
nel [772 (2030) 7t]l = i is produced via 7r 2 (2245); its mass and width are determined better by present 
data than by earlier analyses: M = 2285 ±20 (stat) ±25 (syst) MeV, T = 250 ±20 (stat) ±25 (syst) 
MeV. This is the third channel in which it has been observed. 

Further data on r)37r° from a transversely polarised target would be very valuable. In such 
data, there are interferences between singlet and triplet states. It is likely that such informa- 
tion would allow a complete spin-parity analysis of the data, improving further on the present 
analysis. 
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